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ABSTRACT 


Di  crete  renewal  processes,  until  recent.ly,  have  not  been  applied 
to  the  mathematical  modelling  of  physical  processes.  Analyses  of  such 
renewal  processes  have  proceeded  on  the  basis  of  generating  functions 
but  the  results  are  often  too  complicated  to  be  of  use.  This  paper 
presents  an  alternative  approach  to  discrete  renewal  theory  and  cal¬ 
culates  many  of  the  more  complex  statistics  of  such  processes. 
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I.  INTRODUCTION 


A  renewal  process  is  one  in  which  events  occur  at  time  tj  ac¬ 
cording  to  the  behavior  of  some  underlying  stochastic  mechanism. 

When  the  event  occurs,  the  stochastic  mechanism  is  renewed  (i.e.,  re¬ 
placed  by  a  new  mechanism  identical  to,  but  statistically  independent 
of ,  the  previous  ones .  ) 

The  mathematical  analysis  of  reliability  and  queuing  problems 
has  resulted  in  a  great  deal  of  literature  devoted  to  continuous  time 
renewal  theory  (i.e.,  the  events  can  occur  at  any  time).  Aside  from 
the  volumes  devoted  to  the  applications  of  the  theory,  virtually  all 
books  on  stochastic  processes  devote  at  least  a  section,  if  not  a 
chapter,  to  its  discussion. 

Discrete  time  renewal  processes  (i.e.;  che  events  can  occur  at 
only  fixed  times)  have,  unfortunately,  not  enjoyed  such  a  wide  ex¬ 
posure  either  by  virtue  of  their  apparent  intractability  or  thr  lack 
of  an  obvious  application.  Recently,  however,  applications  ’  have 
arisen  wherein  a  discrete  time  renewal  process  appears  to  be  quite 
useful.  To  be  specific,  the  error  process  in  digital  communication: 
systems  has  been  modelled  by  a  discrete  time  renewal  process.  Anot  lie 
process  being  modelled  by  a  renewal  process,  at  Bell  Telephone  Labo¬ 
ratories  ,  is  the  typing  of  characters  on  teletypewriters.  This  latu 
application  is  quite  important  insofar  as  efficient  communications 
:  y stems  are  being  sought  for  the  remote  programming  of  time-shared 
com put :rs  .  * 


"While  present  typewriter?  are  basically  asynchronous ,  new  Bell  Syste 
models  will  operate  in  a  synchronous  manner. 
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The  discrete  time  renewal  process  would,  no  doubt,  be  quire  use¬ 
ful  in  the  modelling  and  analysis  of  those  systems  in  which  events 
were  fixed  to  occur  at  specific  instants  in  time.  The  growth  of  digi¬ 
tal  technology  has  fostered  the  creation  of  such  systems  in  communi¬ 
cations  and  automata.  It  is  the  purpose  of  this  paper*,  therefore,  to 
investigate  the  statistics  of  such  processes  in  order  to  facilitate 
their  use  in  practical  problems . 

The  statistical  mechanism  governing  the  renewal  events  can  be 
described  in  the  following  manner.  Consider  the  non-negative  random 
variable  X,  which  for  the  purposes  of  discussion,  is  called  the  fail¬ 
ure  time  of  a  component.  This  variable  is  the  length  of  time  between 
renewal  events.  The  distinction  between  the  continuous  time  and  the 
discrete  time  theory  is  made  as  follows : 

(a)  The  random  variable  has  a  continuous  distribution  over  the 
range  (0,®),  its  distribution  being  determined  by  a  proba¬ 
bility  density  function,  f(x).  This  is  the  continuous  case 
and  is  discuss  .  in  great  detail  by  Cox3  and  will  not  be 
examined  here. 

(b)  There  is  a  constant,  T,  such  that  the  only  possible  values 

of  X  are  (T,  2T ,  ...).  The  process  is  determined  by  its 

gap  length  distribution,  p(j),  which  is  the  probability 

that  X  =  jT.  This  latter  ease  is  the  discrete  renewal  proc- 
4  5 

ess  .  Feller  devotes  a  far  from  insignificant  portion  ox 
his  books  to  this  theory,  but  his  work  does'  not  lone  it¬ 
self  to  practical  (e.g.,  actually  calculating  numbers)  ap¬ 
plications  . 

n 

The  approach  taken  by  Feller,  as  well  as  Haight:'",  has'  been 
to  use  the  generating  function  of  X.  This  paper  will  pre¬ 
sent  an  alternative  way  of  looking  at  discrete  time  renewal 
processes  and,  in  particular,  will  present  relatively  simple 
expressions  for  some  of  the  more  complicated  statistics  >1 
such  processes. 

die  work  on  this  paper  was  done  in  the  period  March-June  1  H h . 
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II.  ELEMENTARY  STATISTICS 


The  discrete  renewal  process,  by  definition,  has  gaps  (length  of 
time  between  failures)  whose  lengths  are  independent  and  are  distrib¬ 
uted  according  to  a  common  distribution.  Let  p(j)  be  the  probability 
that  following  a  failure,  the  next  failure  occurs  at  jT.  For  ease  of 
description,  a  trial  is  said  to  be  made  at  every  T.  That  is,  p(j) 

P(0J  "*■  1 1 1 )  where  1  denotes  a  failure,  0  denotes  a  non-failure*  and  0^ 
corresponds  to  i  successive  zeroes.  The  gap  length  distribution  has 
the  following  properties: 


E  p(j)  =  1 

j=i 


CD 


and  that  the  average  failure  rate,  p  ,  is  equal  to  the  reciprocal  of 
the  average  distance  between  failures. 


£  Jpo) 


-1 


(-') 


An  alternative  definition  of  failure  ratt 


Pi  =  lim  cjn^ber  of  xailures  m  (OjJFj 

J-.CE 


'A  non -failure  will  also  be  ca  lied  a  mcces 


Let 


00 

Q(m)  =2  PCd) 

j=m+l 


(3) 


That  is,  Q(m)  is  the  probability  of  there  being  at  least  m  consecutive 
successes  following  a  failure.  Now,  consider  the  event  where  there 
are  m  successes  before  a  failure  (i.e.,  0ml) .  The  probability  of  this 
event  is 


m 

P(0  1)  =  Pl-PlV  p(j)  px  Q(m)  (4) 

3=1 

Note  that  P(0ml)  =  P(1  0m).  This  facilitates  the  calculation  of  the 
probability  of  seme  rather  complicated  events. 

In  many  applications ,  the  probability  of  m  failures  out  of  n 
trials,  P(m,n),  is  important.*  This  is  often  called  the  counting  sta¬ 
tistic.  Elliott7  has  developed  an  equation  for  P(m,n)  for  a  renewal 
process . 


n-m+1 

p(m,n)  Q(j-l)  R(m,n-j+l)  3,imsn  (5) 

3=1 


where  R(m,n)  is  the  probability  that  m-1  failures  occur  in  the  n-1 
trials  following  a  failure.  Therefore,  R(l,n)  =  Q(n-l;  for  n  is  1  and 

n-m+1 

R(m,n)  =  ^  p'j)  R(m-1,  n-j)  2£m£n  ^ 

3=1 

“This  probability  assumes  no  knowledge  of  the  failures  in  the  trials 
preceding  the  n  trial  sequence  under  investigation. 
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It  is  to  be  noted  that  P(m,n)  can  be  expressed  in  terms  oi  and 
p(j),  for  j  <  n .  Several  of  the  above  parameters  and  distributions 

are  analogous  to  those  in  continuous  time  renewal  processes .  For  •  v- 

ample,  p(j)  -  f(x)  and  p,  -  u. 

It  is  at  this  point  that  it  is  worthwhile  to  compare  the  simplic¬ 
ity  of  the  above  expressions  with  those  avaj lable  through  the  use  1 

generating  functions.  Let  g(z)  and  G(z)  be  the  generating  function; 
associated  with  p(j)  and  Q( j ) ,  respectively.  That  is, 

00 

g(z)  =  ^  p( j)  z*  (  ') 

j-i 


and 


=  £  Q(j)  ^  ’ 
j=l 


From  Equation 


d(z) 


i-gyz) 

1-c 


L,<  t  t  it! 


CO 


rv-^rn 


vt  aJtornut  ivc'ly 


H_(2)  =  P,Z 
m  1 


>g(g) 


1-2 


(12) 


Calculation  of  P(m,n),  as  Elliott  puts  it,  is  rather  inconvenient  when 
cor, pared  to  Equation  5.  Generating  functions  do  turn  out  to  be  quite 
useful  for  certain  problems  in  renewal  theory,  but  it  is  noted  that 
the  recurrence  equation  approach  for  P(m,n)  is  ideally  suited  for 
calculations  done  on  a  computer. 

In  addition  to  the  gap  length  distribution,  the  "autocorrelation" 
of  the  failures  is  often  quite  useful.  The  autocorrelation,  a(j),  is 
defined  as  the  probability  that,  following  a  failure,  a  failure  occurs 
at  X  ~  jT.  For  a  renewal  process,  however,  the  autocorrelation  can 
be  derived  from  the  gap  length  distribution. 

j-1 

a(j)  =  P(j)  +  ^  p(s)  a  (j-s)  for  j>l  (13) 

s=l 


where  a(l)  =  p(l)  and  a(0)  =  1.  Therefore  a  renewal  process  can  be 

4 

described  by  p(i),  a(-j)  or  P(m,n).  Felle^  has  also  proven  that 


lim 

j-  00 


a(j)  =  P1- 


A  renewal  process  c  interest  is  one  in  which 

a(j)  =  a8“  +  Px  j  ^  1  8  <  1  .  (14) 

8 

This  is  the  autocorrelation  Gilbert  obtains  in  his  Markov  model.  It 
happens  that  this  simple  model  generates  a  renewal  process  while  the 
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Berkovits  model,  which  has  the  same  autocorrelation ,  is  not  a  renewal 
process.*  For  a  renewal  process  with  this  autocorrelation,  and  if 
p^  «  0,0  for  j  £  N ,  then 


a(j)  «  oBJ 


From  Equation  13 


p(j)  =  o8J(l-o)j_1 


Then  if 


oB 


1-8  (1-0 0 


1, 


Q(m) 


_  tf"*1  (1-g 

TT-ba 


1-a), 


(if) 


(16) 


(17) 


and 


,  b1]  am  8n(i^)n 

R(m,n)  -  \p_xj  Tl-3(1^TT 


-m 


Using  Equation  5 

P(m:n)  -  p. 


0 


gmtl  8n+1  (l<On-m 

n-8(l-c012 


1  5  m  5  n 


(18) 


(19) 


At  this  point,  it  is  apparent  that  a  Bernoulli  process  (i.e., 
one  in  which  each  trial  is  independent  of  the  others)  is  a  renewal 
process  with 

T(j)  =  P  » 

P(j)  =  P(l-P)^"3'  >  (;1') 

/n\ 

\  (  1  m  f-.  m-m 

P(m,n)  -\raj  p  (1-p) 

"The  loot  that  processes  have  the  same  autocorrelations  does  not  imply 
that  tr.ey  have  the  same  gap  length  distributions  unless  thev  are  both 
rt  r.ewal  processes  . 
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III.  COMPLEX  STATISTICS 


A.  COMPOUND  COUNTING  DISTRIBUTION 

The  fact  that  there  may  be  a  correlation  between  failures  raises 
the  importance  of  the  compound  counting  statistic,  PCrn^rn^ ,  *  .  •  5ny,  ,N) . 
This  probability  corresponds  to  the  following:  Consider  a  sequence 
of  mN  trials,  as  shown  in  Fig.  1,  which  are  divided  into  N  subsequences 
of  n  trials  each.  In  the  1  subsequence  there  are  rru  failures 

J  <  m.  s  n 
1 

The  first  failure  in  the  i^  subsequence  is  preceded  by  -  1  suc¬ 
cesses  anc  the  last  failure  is  followed  by  a .  -  l.  successes.  The 

l  l 

failures  are  clustered  in  an  interval  of  length  n  -  +  1.  Thi^  con¬ 

tinues  until  the  last  n  trials  where  ••  1  failures  occur  in  the  last 
m  -  i  trials  though  not  necessarily  ending  in  a  failure.  The  black 
bars  show  the  beginning  and  end  of  an  rru  failure  cluster  in  each  n 
trial  subsequence.  To  examine  this  case,  first  let  S(m,n)  be  the  prob¬ 
ability  that  (m  -  1)  failures  occur  in  the  (n  -  1)  trials  following  a 

Q  t  c  t 

failure,  and  that  the  (m  -  1)—  failure  is  at  the  (n  -  1)—  trial. 
S(m,n)  obeys  the  following  recurrence  relation 


n-rn+1 

S(m,n)  =  Y]  S(m-1,  n-j)  p(j)  (.'1) 


for  2  <  in  <  n  where  S(l,l)  =  0  and  S(l,n)  =  0  for  n  >  I.  R(m,n)  in¬ 
troduced  in  Equation  5  follows  directly  from  S(m,n)  since 


B 


(22) 


n-m+1  oo 

R(m,n)  =  ^  ^2  S(m,n- j+1)  p(i) 

j=l  i=j 


FIGURE  1.  Failure  Sequence  for  Complex  Counting  Statistics 


The  general  expression  for  the  compound  counting  statistic  is 


P(m1.m2....,mN,N)  -  Px  £ 


Q  (^-D 


N-l 


TTS  ( m  . ,  n-a  +1 )  p ( a -I  .  +t  +±+l )| 


1=1 


(23) 


R(mN  ,n-<tN+l) 


where  nrn  >  0.  The  first  (N-l)  pairs  of  summations  over  a.  and  l 
summed  as  follows: 


are 


n-m.+l  n-m.+l 

£ 

t  .=1  a  ,=-t  . 

1  2  1 


for  1  s  j  sN-1 


n-mN+l 


and  tne  last  summation  is 


buppose  m 


0  (i.e.,  a  group  of  trials  has  no  failures).  This 


can  happen  in  three  ways: 


(i) 


(ii) 


(iii) 


If  the  n  trial  block  is  the  first  block,  omit  S(m1 ,  n-a^+1) 
and  pCa^^+^+l)  and  the  associated  summations  and  change 
Qtt^l)  to  Q(t1+n-l). 

If  the  n  trial  block  is  the  last  block,  then  omit  the  term 
R(rrL.,n~t..+l)  and  sum  the  last  term  p(a..  , -llT  n+.Lt+l)  over 
the  limits  m  <  <  00 . 

If  the  n  trial  block  is  internal,  omit  the  respective 

S(n.,  n-a  .+1)  p(a  .-l  .+t  ■ , , +1)  term  as  well  as  the  corre- 
3  J  J  3  3+1 

s ponding  summations  over  a_.  and  i  ^  and  change  p(a_.  ^ 


+1  .+1)  to  p(a . 

3  J- 


.  -i .  ..  +1 .+  n  +  1) 

■1  j-l  j 


The  above  process  of 


removing  terms  can  be  cont-i'.nued  up  to  the  point  where  only 
0  because  Equation  23  assumes  at  least  one  renewal 


one  m_ 
event . 


For  all  m.  =  0, 


P(0  , . .  .,0,N)  =  p,  X)  QCj) 

j=mN 


B.  AGING 

An  interesting  statistic  often  calculated  for  continuous  renewal 
processes  is  the  age-specific  failure  rate.  This  is  defined  as  the 
probability  of  a  failure  between  x  and  x+Ax  given  it  has  not  failed  up 
to  x.  For  the  discrete  case  it  is  defined  as  he  probability  of  a 
1  allure  on  the  j —  trial  given  no  failures  up  to,  and  including,  the 

c,  +- 

( j -1)- —  trial  with  no  knowledge  of  the  trials  preceding  these  j  trials. 
Denoting  it  by  cp(j),  one  gets 


cp(j) 


QC.i-D 

£  Q(m) 
m=j-l 


(21) 


% 
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If  cp(j)  increases  with  j  it  is  said  to  have  positive  aging  and 
becomes  more  likely  to  fail.  Some  processes  have  a  ?(  )  that  decreases 
with  j.  There  is  however  a  process  with  no  aging.  Tha :  is,  let 


cp(j  )  =  X 


or 


(25) 


Q(j-l)  =  X  £  Q(m)  • 
m=j-l 


Noting  that  Q(0)  =  1  and  53  Q(nO  = 

m=0 


—  vields 

Pi  ' 


X  = 


P1 


and  (26) 

p(j)  =  Pl  (1-p^-'"1  , 

which  is  obviously  the  Bernoulli  process. 


C.  BURSTS  OF  FAILURES 

A  simple  extension  f  the  previous  discussion  is  to  consider 
bursts  of  failures.  A  burst  of  length  m  is  defined  to  be  m  consecu¬ 
tive  trials  beginning  and  ending  a  failure.  This  does  not  require 
all  m  trials  to  be  failures,  but  unly  the  first  and  last .  When  con¬ 
sidering  n  consecutive  trials,  the  probability  of  it  containing  a 
burst  of  length  m  is  denoted  by  B(m,n).  Bursts  of  faili  rer  are  of 
considerable  interest  in  many  physical  processes  which  may  be  modelled. 
For  example,  errors  in  digital  communications  are  not  ed  to  be  clu:  terc 
(i.e.,  come  in  bursts)  and  error  correcting  codes  are  constructed  t 


deal  with  these  bursts.  It  is  also  noted  that  the  hittinq  of  tv; 
writer  keys,  which  in  this  case  corresponds  to  failure,  occurs  in 
The  examination  of  bursts  of  failures  is  therefore  important. 


bur: 
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For  the  process  described  previously 


B(m,n)  =  p  [a(m-l)j  [p(n-m+l)  +  p(n-m+2)+. . . ] 

+  [l-p(l)l  p1  [a(m-l)]  [p(n-ra)+p(n-m+l)+. ..]  (27) 

+. . .+[l-p(l)-p(2)-. . .-p(n-m)]  p1  [a(m-l)] 

The  first  term  in  the  above  expression  corresponds  to  the  first  fail¬ 
ure  of  the  burst  in  the  first  of  the  n  trials;  the  second  term  to 

where  the  first  failure  is  at  the  second  trial;  and  the  last  term  is 

C  t  — 

where  the  first  failure  of  the  burst  is  on  the  (n-m+1)—  trial.  This 
can  be  compactly  written  as 

n-m 

B(m,o)  =  p,a(m-l)£  Q( j )Q(n-m- j)  (28) 

j~0 


D .  INTERLEAVING 

A  technique  often  used  to  overcome  the  effects  of  bursts  is  to 
interleave  the  process.  For  example,  in  the  case  of  digital  trans¬ 
mission,  time  division  multiplexing  M  data  streams  such  that  two  dig¬ 
its  originally  adjacent  in  any  one  stream  are  transmitted  M  digits  a- 
part  spreads  a  burst  of  errors  over  the  M  streams  and  makes  the  errors 
look  almost  as  if  they  were  randomly  distributed  in  any  one  data 
stream.  For  the  remote  programming  application,  M  parallel  computer 
input  ports  would  be  able  to  digest  a  burst  of  telety pewriter  charac¬ 
ters  if  each  port  accepted  every  character.  The  pro’'lem  is  now 
to  analyze  the  process  obtained  by  examining  every  M  '  trial  of  a  dis¬ 
crete  renewal  process. 

The  autocorrelation  of  the  events  in  the  interleaved  process, 
a(j,M),  is  obviously  a(jM).  The  gap  length  distribution  is  now  given 
by 
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(29'' 


P(j,M)  =  a(jM) 


j-1 

X)  a(sM)  P(j-s,M)  j  >  1 
s=l 


As  M  gets  large, 


lim  a( jM)  =  p.  (. 


Therefore 


lim  p(j,M>  -  p1(l-p1)3"1  (31) 

M-*°°  *  'L 


and  the  process  becomes  a  Bernoulli  process. 

E.  FREEZE- O' !T  PROBLEMS 

The  freeze-out  problem  arises  when  a  physical  device  requires  a 
specific  time  interval  to  digest  recorded  data  and  cannot  accept  ad¬ 
ditional  data  during  that  time  interval.  For  example,  an  event  occur 
randomly  in  time.  When  it  occurs  a  "counter"  must  record  some  appro¬ 
priate  data.  This  recording  period  lasts,  say,  n  seconds  (or  trials) 
and  if  another  event  occurs  during  these  n  seconds  it  goes  unrecorded 
The  probability  of  such  an  event  for  a  discrete  renewal  process  being 
unrecorded  is  not  examined  directly  but  rather  two  different  measures 
are  used:  the  mean  time  to  ar.  unrecorded  event  and  the  probability 
ol  at  least  one  unrecorded  event  in  N  trials. 

1 .  Mean  Time  to  Unrecorded  Event 

If  the  probability  that  the  first  unrecorded  event  occurs  on  the 

t.  h  — - 

j —  trial  is  denoted  by  P_. ,  the  the  mean  time  to  failure,  T,  is  nave 

by 
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1 

1 

J 


The  guard  space  time  to  prevent  a  freeze -out  is  n-1  trials  such  that 
if  two  events  occur  within  any  n  trials ,  then  the  second  is  not  re- 
cor^ed.  In  the  analysis  that  follows,  the  signal  flow  graph  techniques 
of  Sittler10  and  Huggins11  will  be  used.  The  state  diagram  of  Fig.  2 
is  obtained  where : 

i  is  the  state  of  having  at  least  i  successes  (in  this  case 
trials  in  which  events  do  not  occur)  before  the  first  failure  (event); 

i  is  the  state  of  having  at  least  i  successes  since  the  last 
failure ;  and 

is  the  state  of  having  a  failure  but  not  an  unrecorded  event. 

Then  letting  q(a,b)  be  the  probability  of  going  to  state  b  from  state 

a, 


qd.TJT) 


1-P,£n  Q(m) 

- FT - 

1*PiE  Q(m) 

m=0 


(33) 


P-,  Q(i) 

qd.E^  -  - ^ -  (34) 

1_Pi  £  Q(nO 
m=0 


Q(U12 

QU) 


(  3‘  ) 


q(i,Ex) 


1) 


for  i  -  n-1 


(  ) 
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for  i  £  n-2  (37) 


q(_i,  unrecorded  event)  = 


The  flow  graph  can  be  reduced  to  that  of  Fig.  3.  The  "transmission 
from  start  to  freeze-out"  is  P(Z). 


Pxz 

00  m 

£  Q(m)  Zm 

-m=0 

r  n- 1  .q 

E  P(3) 

-i=l 

P(Z) 

« 

1 

-E  p(j) 

j=n 

If  P(Z)  were  expanded  as 

P(Z)  - 

P^+PXZ+P2Z2 

+  . . .  +P  -Z-*  i  •  •  • 

3 

the  coefficient  of  ZJ  is  the  probability  of  an  unrecorded  event  on  the 


dC 

Q  (m)  -  Z  p(j )  I  -  ;  p{l  )  ♦  .  . .  ♦  p(m) ; 

j  m  +  ! 

FIGURE  3.  Reduced  Stote  Flow  Graph  for  a  Freeze-Ouf  of  i  Counter 
Requiring  a  Guard  Space  of  n-1 


From  the  definition  of  T, 


f  =  Mil 

dZ 


Z=1 


(40) 


and  therefore 


T 


rWn-iy] 


+  Pi 


m  Q(m) 


(41) 


2 .  Freeze-Out  in  N  Trials 

The  preceding  analysis  evaluated  the  mean  time  to  freeze-out, 
but  this  measure  may  have  little  significance  to  someone  designing  an 
experiment  to  record  data.  As  an  alternative,  P(N),  the  probability 
of  at  least  one  freeze-out  in  N  trials  will  be  examined. 


where 


N 

P(N)  =  D  P. 
j=l  3 


(42) 


(43) 


The  complicated  nature  of  P(Z)  precludes  this  approach. 


Another  technique 
the  only  way  to  have  a 
failure  in  the  first  j 


is  to  calculate  P.  directly.  For  j  £  n  +  1, 

1  th 

freeze-out  on  the  j —  trial  is  to  have  only  one 
-1  trials  .  Then 
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QOO  Px  p(j-k-l) 


(44) 


k=0 


who  is.  the  above  expression  has  k  successful  trials  before  the  first 

A.  Vs 

failure.  If  n  +  1  <  j  s  2n  +  1,  a  freeze-out  on  the  j —  trial  can 
occur  because  of  a  failure  in  trials  j  -  n  +  1  to  j  -  1,  or  by  two 
previous  failures.  The  probabilities  are 


Q(k)  P1 


p(.i-k-l) 


(45) 


j-n-2  3-k-2 

and  ]T  22  Q(k)  P(r>  P(j-k-r-l)  ,  respectively. 

k=0  r=n 


Continuing 
that  if  in 
l  .  a r>  prior 


in  this  vein  becomes  extremely  difficult  when  one  realizes 
+  1  <  j  s  (i  +  1)  n  +  1  then  it  is  possible  to  have  m  fail- 
to  the  failure  causing  the  freeze-out  where  m  =  1,  ..., 


i  -»  1. 


An  alternative  approach  is  to  work  backwards.  It  was  shown  in 
Equation  4  that  reversing  a  pattern  of  failure  does  not  alter  the 
probabilities  involved.  Then  for  a  freeze-out  on  the  j — -  trial,  there 
can  be  at  most  n  -  2  successes  preceding  it  after  the  previous  failure. 
Preceding  that  failure  there  is  at  least  n  -  1  successes  and  so  on. 

The  probability  of  a  freeze-out  on  the  j — -  trial  is  now  compactly  ex¬ 


pressed  as 


P. 

J 


n-1  j-(m-l)n 

E  E 

rQ=l  r1=n 


j-(m-2)n-r 

2 


;-n-r, ...-r  _ 

1  m-2 


r  =n 
m-1 


hp(r0} 


P(ri>. .  .ptr^)  Q(  j-r0-rr  . .  -r  -1) 
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(46) 


where  rn  is  defined  above.  The  only  exception  to  the  above  occurs 
when  m  =  i  +  1  and  j  =  in  +  a  where  1  <  a  <  n  in  which  case  Equation 
46  becomes 


a-1 

Pj  =  Px  I’pCn)]1  p(X)  Q(j-ro-l-in)  (47) 

r  =1 

o 

Here  again,  as  in  the  previous  statistics,  a  complicated  probability 
is  obtained  in  a  form  directly  amenable  to  computer  calculation  where¬ 
in  the  only  data  needed  is  p^  and  p(j)  for  j  s  N  -  2. 
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